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OPERATING  CHARACTERISTICS  FOR  MAXIMUM  LIKELIHOOD 
DETECTION  OF  SIGNAL  IN  GAUSSIAN  NOISE 
OF  UNKNOWN  LEVEL 

H.  PHASE-INCOHERENT  SIGNALS  OF  UNKNOWN  LEVEL 


INTRODUCTION 

In  the  first  report  of  this  series  on  the  operating  characteristics  for 
detecting  unknown -level  signals  in  Gaussian  noise  (reference  1),  maximum 
likelihood  (ML)  detection  of  signals  of  known  form  in  additive  Gaussian  noise  of 
unknown  level  was  analyzed.  For  the  case  of  coherent  signals  corrupted  by 
noise  of  unknown  level,  it  was  shown  that  the  ML  detector  is  optimum  in  the 
sense  that,  of  all  processors  yielding  a  specified  false-alarm  probability  Pp- 
without  knowledge  of  absolute  levels,  the  ML  detector  is  uniformly  most 
powerful. 

In  this  report,  the  previous  study  is  continued  to  cover  two  cases  of  ML 
detection  of  nonfluctuating,  phase-incoherent  signals  in  noise  of  unknown  level. 
Again,  as  in  the  case  of  the  previous  study,  the  ML  detectors  will  be  shown  to 
possess  the  important  property  of  maintaining  a  specified  Pp1  irrespective  of 
the  actual  noise  level.  However  a  major  difference  exists  between  the  present 
and  previous  analyses.  Whereas,  for  the  coherent  signals,  the  ML  detector  is 
optimum  for  the  class  of  detectors  invariant  to  level  scaling,  no  such  conclusion 
is  possible  for  the  phase-incoherent  signals.  Because  of  this  inability  to  resolve 
the  question  of  performance  quality,  we  shall  consider,  but  not  thoroughly 
analyze,  several  reasonable,  albeit,  ad  hoc  processors.  That  is,  we  shall 
discuss  and  compare  the  performances  of  ML  processors  with  other,  heuris- 
tically  obtained,  processors.  Relevant  past  work  on  this  problem  is  furnished 
by  reference  2  for  the  special  case  of  a  single  signal  sample. 

In  the  analysis  that  follows,  it  is  supposed  that 

1.  sufficient  knowledge  of  the  signals  is  available  to  allow  filtering  the  in¬ 
put  to  the  receiving  system,  thereby  removing  noise  outside  the  signal  band,  and 

2.  the  sampled  outputs  of  the  system  filters  are  statistically  independent. 

In  order  to  restructure  the  originally  nonparametric  problem  into  a  parametric 
form,  it  is  further  assumed  that  the  noise  is  Gaussian  and  stationary. 
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In  determining  the  ML  detectors,  the  values  of  the  unknown  parameters 
are  estimated  by  separately  maximizing  the  probability  density  functions  (PDFs) 
for  the  signal-present  hypothesis  and  the  signal-absent  hypothesis  Hq.  It 
is  reasonable  to  expect  that  the  noise-level-estimation  procedure  is  improved 
if  extra  samples  consisting  of  noise-alone  are  available.  This  possibility 
which  is  of  practical  importance,  is  considered  here.  Indeed,  in  one  case,  the 
ML  detector  has  meaning  only  if  extra  noise-alone  samples  are  available. 

The  likelihood  ratio  (LR),  formed  using  the  ML  estimates  in  the  PDFs, 
results  in  a  generalized  likelihood  ratio  (GLR)  test.  In  this  report,  the  GLR 
procedure  is  used  for  two  cases  involving  different  assumptions  about  the  signal. 
In  the  first  case,  the  signal  amplitudes,  when  signal  is  present,  are  assumed  dif¬ 
ferent  and  unknown  from  sample  to  sample.  This  will  be  referred  to  as  the  unequal 
amplitude  (UA)  case.  In  the  second  case,  the  signal  amplitudes,  although  unknown, 
are  assumed  to  be  the  same  for  all  samples.  This  will  be  referred  to  as  the  equal 
amplitude  (EA)  case.  In  both  cases,  the  signal  amplitudes  are  nonfluctuating. 

The  performance  characteristics  of  the  ML  detector  for  the  UA  case  can 
be  obtained  analytically,  regardless  of  the  actual  distribution  of  signal  ampli¬ 
tudes,  and  extensive  results  are  presented.  For  the  EA  case,  these  perform¬ 
ance  characteristics  are  not  amenable  to  closed-form  solution,  however,  and 
receiver-performance  simulation  has  been  conducted. 

In  the  following  sections,  the  problem  is  stated  in  detail,  and  ML  detectors 
for  the  UA  and  EA  cases  are  derived  and  discussed.  The  performance  charac¬ 
teristics  of  the  UA  processor  are  derived,  and  a  description  of  the  simulation 
procedure  for  the  EA  processor  is  given.  To  ascertain  the  usefulness  of  the 
ML  detectors,  the  detection  performances  of  several  ad  hoc  processors,  whose 
Pp  is  invariant  to  noise-level  scaling,  are  also  briefly  considered. 


PROBLEM  STATEMENT 


It  is  supposed  that  a  detection  system  performs  M  observations  of  a 
received  signal-plus-noise  or  noise-alone  process  having  waveforms 


Re{Vk(t)  exp(i2xfkt)[ 


1  <  k  <  M 


(1) 


where 


1  <  k  <  M  . 
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In  (2),  S^t)  and  N^(t)  are  the  complex  envelopes  of  the  received  signal  and 
Gaussian  noise  waveforms,  respectively,  for  the  k-th  observation,  and,  in  (1), 
f^  is  the  center  frequency  used  in  the  k-th  observation.  We  note  that  the  M 
observations  referred  to  above  may  denote,  for  example,  successive  observa¬ 
tion  intervals  in  time,  simultaneous  observations  either  in  different  frequency 
bands  or  in  different  spatial  locations,  or,  indeed,  combinations  of  these.  These 
M  observations  are  called  the  potential- signal  observations. 

Further,  it  is  assumed  that  N  additional  observations  are  available  in 
which  it  is  known  that  only  noise  is  present.  For  these  latter  observations,  the 
received  waveforms  are 


Re  jv.(t)  exp(i2?rf  t)}  ,  1  <  j  <  N  , 

where 

V.(t)=N.(t),  l<j<N, 


(3) 

(4) 


for  both  and  H0 . 

Following  standard  procedures  (reference  3),  receiver  processing  con¬ 
sists  of  computing  the  complex  quantities 

Jdt  Fk(t)  Vk(t)  =  zk  =  xk  +  iyk,  1  <  k  <  M  ,  (5) 

and 

fdt  F.(t)  V.(t)  5  2.  =  x.  +  iy.  ,  1  <  j  <  N  .  (6) 

(Integration  is  over  the  range  of  nonzero  integrand. )  {F^(t)}  and  {Fj  (t)}  are 
filters  that  need  not  be  optimum  (matched)  in  any  sense.  This  generality  in 
filter  specification  allows  for  mismatch  in  frequency  or  in  time,  or  both, 
between  transmission  source  and  receiver.  Because  {^(t)}  and  (Rfj(t)}  are 
Gaussian  processes,  {zk}  and  {2j}  are  complex  Gaussian  random  variables. 
They  are  proportional  to  samples  of  the  complex  envelopes  at  the  filter  outputs. 

At  this  point,  we  restrict  attention  to  the  case  where  the  M  +  N  filter 
outputs  are  statistically  independent  of  each  other  and  of  equal  (but  unknown) 
variance  (see  appendix  A).  The  above  restriction,  for  example,  would  be 
pertinent  when  the  receiver  filters  are  disjoint  in  time  or  frequency  and  of 
equal  gain,  and  when  the  noise  spectrum  is  substantially  flat  over  the  entire 
detection  system  filter  bank. 
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To  summarize  the  discussion  thus  far,  we  are  considering  an  observation 
vector  R  of  independeht  Gaussian  random  variables: 

rT  yi‘  ^mVi  "  Vn]'  <7> 

On  the  basis  of  R,  we  are  to  decide  between  hypotheses  (signal  present 
in  and  H0  (signal  absent  in  {x^.y^}).  The  variance  associated 

with  R  is  imknown;  nevertheless  we  wish  to  consider  processing  of  R  that 
will  ensure  a  specified  Pp  irrespective  of  the  magnitude  of  As  already 
noted,  the  ML  technique  is  used  to  estimate  all  unknown  parameters,  and  the 
GLR  (references  3,4)  is  used  to  determine  a  test  statistic.  This  procedure  is 
discussed  in  the  next  section. 


ML  ESTIMATION  PROCEDURE  AND  ML  DETECTOR 


Under  hypothesis  H0 ,  the  PDF  of  observation  R ,  conditioned  on  an 
assumed  value  of  noise  variance  ,  is  (appendix  A) 


>>2)=n 


k=l 


exp 


2  2^. 

Vyk 


2  ira 


2  a 


N 


n 

jj  j=i 


exp 


r.ii 


2ir(j 


2  <T 


(8) 


The  ML  estimate,  denoted  by  carats,  of  oq  imder  hypothesis  H0  is  that  value 
of  <r^  which  maximizes  (8),  namely, 

M  N 


*2  k=l 

2  9  - - 


j=l  k=l 


t 

j=l 


M  +  N 


M  +  N 


(9) 


which  is  proportional  to  the  total  sample  power  over  all  M  +  N  samples. 
Substitution  of  (9)  in  (8)  yields 


M 


max  p 
2  o 

a 


(rU2)  = 


e  r  ‘ 


,  ,2  ^  H2\-<M  +  N> 
ki  +  z;  H 

k=l  J=1 

M  +  N 


z 


(10) 
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Under  hypothesis  Hi ,  the  PDF  of  R  ,  conditioned  on  assumed  values  of 
signal  means  and  ,  is  (appendix  A) 


v  M 

L2  a  ii/)  —  T"T 

1 

'  cos  *k)2 + (yk  -  \ sin  •‘‘tS  ] 

k=l 

r  exp 

2xa 

2<r2  J 

(11) 


N 

n 

j=i 


~2 


exp 


2ito 


*s 


2a 


where,  denoting  ensemble  averages  on  the  noise  by  overbars,  the  means  are 
(from  (A -4)) 


and 


zk  =  xktiyk2Akexp(i*k)1  1sk±M’ 


<M[V"  *Mf  • 


(12) 


(13) 


The  ML  estimate  a“  obtained  under  hypothesis  is  (from  (11)) 


M 


N 


2S2  =  — 


ElvW^ZN 

j=l 


M  +  N 


Substituting  (14)  in  (11),  we  obtain 


(14) 


max  p1(R|<r2,A,l//)  = 


M  .2  Ai  ,2\'(M  +  N) 

El\ - \ exp(i*k)i  +Ehl N 


j=i 


M  +  N 


(15) 


This  function,  which  pertains  to  hypothesis  Hi ,  can  be  further  maximized  by 
choice  of  assumed  signal  amplitudes  A  and  phases  \f/.  To  proceed  with  this 
program,  we  distinguish  two  cases  pertaining  to  different  assumptions  about 
the  signal  amplitudes  {A^}.  The  first  case,  UA,  will  apply  to  the  situation 
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where  each  {Ak}  is  allowed  to  be  different;  the  second  case,  EA,  will  apply  to 
the  situation  where  each  {Ak}  is  assumed  to  be  the  same  unknown  value;  that 
is,  Ajc  =  A,  l<k<M. 


For  the  UA  case,  the  ML  estimates  of  signal  amplitudes  and  phases  are, 
upon  use  of  (15), 

=  |zfe |  »  \  =  arS(zk)  »  1  <  k  <  M  .  (16) 

Substitution  of  the  ML  estimates,  (16),  in  (15)  gives 


max  p 

*2,a,<|/ 


(r|«2,a,\|/)  = 


A  GLR  test, 


max  p1(rU2,A,v(/) 

*  ,A„lJ/  _ 

max  p  (RU  ) 

2  o' 

a 


>  threshold , 
H 

o 


(17) 


(18) 


where  hypothesis  is  declared  true  if  the  upper  inequality  is  satisfied,  and 
hypothesis  H0  is  declared  true  if  the  lower  inequality  is  satisfied,  can  now  be 
constructed.  For  the  UA  case,  GLR  test  (18)  is  expressible  as 


M  >  1,  N  >  1  , 


(19) 


where  r(>0)  is  a  constant.  The  test  states  that  the  sums  of  squared-envelope 
samples  at  the  filter  outputs  for  potential-signal  samples  and  noise-alone  sam¬ 
ples  be  compared. 

Test  (19),  which  when  implemented  is  the  UA  processor,  depends  explicitly 
on  the  availability  of  noise -alone  samples;  that  is,  if  N  =  0  ,  a  test  different 
from  the  one  above  must  be  formulated.  It  is  noted  that  implementing  the  UA 
processor,  (19),  results  in  a  constant  false  alarm  rate  (CFAR)  receiver,  since 
both  sides  of  the  test  can  be  scaled  arbitrarily  for  noise-alone  without  changing 
the  result  of  the  test. 
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The  EA  case  can  be  handled  similarly.  The  ML  estimates  of  the  signal 
amplitude  and  phases,  which  are  available  from  (15),  follow: 

M 

s=¥Lk|’  <2°> 

k=l 


which  is  the  average  envelope  for  the  potential-signal  samples,  and 

a 

=  arg(z^)  ,  1  <  k  <  M  . 

Substituting  ML  estimates  (20)  and  (21)  in  (15),  we  obtain 


(21) 


M  .  M 

Em "ff(£h 


12  I 

a  ,  A,\J/)  =  leir- 


^k=l 


The  GLR  in  this  case  then  becomes 
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M  +  N 


-(M  +  N) 


GLR  = 
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k=l 


N 

♦Eh 

j=l 
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M  +  N 


‘Jl 
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M 


,2  N 


Ekl  24(EhJ)  +E 

\  k=l  k=l  j=l 

and  the  GLR  test  can  be  written  in  the  form 
M  \2  /  M  N 


I  —  I  2 
z. 


Ekl  (El1 


k=l 


k=l 


+E 

j=i 


~  ,2' 

Zi 


M  >  1,  N  >  0 


(22) 


(23) 


(24) 


where  7  is  a  threshold  (scale  factor).  The  left  side  of  (24)  is  the  square  of 
the  sum  of  envelopes  at  the  potential-signal  filter  outputs.  The  right  side  is 
proportional  to  the  sum  of  squared  envelopes  of  all  the  filter  outputs. 


Test  (24),  which  when  implemented  is  the  EA  processor,  is  a  CFAR  test. 
It  becomes  meaningless  for  the  situation  M  =  1,  N  =  0,  and,  of  course,  no 
realization  of  a  CFAR  processor  is  possible  in  this  instance.  The  test  is 
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similar  in  form  to  the  GLR  test  obtained  for  phase-coherent  signals  in 
Gaussian  noise  of  unknown  level  (reference  1).  The  difference  between  the 
tests  for  the  phase-coherent  and  the  phase-incoherent  cases  is  that,  in  the 
former  test,  phase  information  is  preserved  in  the  samples,  while  in  the  latter 
test,  (24),  only  the  sample  envelopes  are  used. 

We  may  note  one  further  point:  For  the  situation  M  =  1  (N  >  1) ,  the  UA 
processor,  (19),  and  the  E A  processor,  (24),  treat  the  available  observations 
identically.  This  observation  is  important  because  not  only  is  the  M  =  1  situa¬ 
tion  of  considerable  practical  interest,  but  also,  as  will  be  shown  in  the  next 
section,  performance  for  M  =  1  for  the  EA  processor  is  thereby  amenable  to 
a  complete  analytical  treatment.  (The  UA  processor  analysis  is  tractable  for 
all  values  of  M  and  N . ) 


EVALUATION  OF  UA  AND  EA  PROCESSOR  PERFORMANCES 


The  evaluation  of  the  performance  of  the  GLR  test,  (19),  for  the  UA  case, 
although  tractable,  is  lengthy,  and  the  derivation  is  carried  through  in  detail  in 
appendix  B.  One  form  for  the  probability  of  detection  Pq  obtained  in  appen¬ 
dix  B  is 


p  =l-( 

f  a  > 

D  \ 

where 

1*1 

M 


expi- 


-l+k\  ,  /  d 


— ) 
1  +  a/ 


,(25) 


total  signal -to -noise  ratio  (power)  parameter  d^/2  appearing  above  is  given  by 


M  M 

K-XKKEi 


Jdt  Fk(t)  Sk(t)|2/2 


k=l 


fel  f  J/dt  du  Fx(t)  F*(u)  Ru(t  -  u) 


(26) 


where  Sjc(t)  is  the  actual  received  signal  complex  envelope  on  the  k-th  obser¬ 
vation,  and 


Nk(t)  N*(u)  =Rkfl(t  -  u)  ,  1  <  k,4  <  M  ,  (27) 

is  the  actual  complex-envelope  noise  crosscorrelation  on  the  potential -signal 
samples.  In  (26),  F^(t)  is  the  complex  envelope  of  the  receiver  filter  em¬ 
ployed  on  the  k-th  potential-signal  branch. 
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For  the  special  case  of  white  noise  and  equal-gain*  matched  filters, 


(28) 


where  Ek  is  the  received  signal  energy  on  the  k-th  observation  and  N(j  is  the 
double-sided  noise  power  density  spectral  level.  In  this  instance,  there  follows 


(29) 


and 


(30) 


where  N0  is  the  single-sided  noise  power  density  spectral  level.  The  quantity 
dj^2  can  be  interpreted  as  the  output  signal-to-noise  ratio  of  the  filter  on  the 
k-th  potential-signal  branch;  it  is  proportional  to  the  received  signal  power  on 


ence  in  the  required  decibel  level  between  two  processors  (for  example,  different 
M  and/or  N)  measures  directly  the  difference  in  input  signal-to-noise  ratio  in 
decibels  required  for  the  two  processors  to  realize  the  same  Pp  and  Pj^ . 

From  (30),  we  note  that  the  exact  way  that  the  total  received  signal  energy 
Ep  is  actually  fractionalized  into  the  observation  branches  is  not  important; 
only  Et  affects  UA  processor  performance.  Thus  the  UA  processor  is  robust 
with  respect  to  actual  signal  energy  fractionalization;  it  represents  an  attain¬ 
able  level  of  performance  regardless  of  the  exact  conditions. 


♦The  normalization  in  (28)  makes  the  filter  gains,  as  measured  by 
/  dt  |Fk(t)|2  ,  independent  of  k  and  makes  <r2  independent  of  k,  as  required 


by  (A-13). 
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When  dp  is  set  equal  to  zero  in  formula  (25)  for  Pjj  ,  we  obtain 

/  VM  N-l  /M  -  1  +  k\  , 

’-•-(*) £f  .V-  • 


(31) 


lk=0 

An  alternative  form  for  Pp  is  obtained  in  appendix  B  and  is  given  by 

pD  =  1'-'rf^  “»(-t  rb)Z  <1+“>'k  !Fi(-k;1;-T  rh) 


k=0 

(1  +«fN  exp(-^  f  1+k)(rfs)lt|1Fi(1-N;k+1;-T;  if;)  • 

k*i'  k  /  \ 

(32) 


(33) 


For  dx  =  0 ,  we  obtain  from  (32) 

m^I/m  - 1  +  k 

PF  =  <1+“»  E(  „ 
k=0  '  ^ 

Comparison  of  the  alternative  expressions  for  Pp  indicates  that,  computa¬ 
tionally,  (31)  is  more  efficient  for  M  >  N ,  whereas  (33)  is  better  suited  for 
N  >  M .  Also  (33)  possesses  only  positive  terms. 


)(*)’ 


For  the  UA  processor,  as  N-*  «»  ,  the  expressions  (25),  (31),  (32),  and 
(33)  for  Pjj  and  Pp  tend  to  the  corresponding  forms  obtained  for  envelope- 
squared  detection  with  known-noise  level  (reference  3),  as  indeed  they  must, 
because  ML  estimators  are  consistent  (see  appendix  B). 


Let  us  consider  now  the  EA  processor.  It  has  already  been  noted  that  the 
EA  and  UA  processors  are  identical  for  the  situation  M  -  1  and  any  N  >  1.  Thus 
the  Pjj  and  Pp  formulas,  assembled  above,  also  are  applicable  to  the  EA  proc¬ 
essor  performance  for  M  =  1.  One  other  analytical  result  for  the  EA  processor 
has  been  obtained:  In  appendix  C,  it  is  shown  that  for  M  =  2  and  arbitrary  N 


r(f)  T(N  +  2) 


N  +  l/2 


,  3\  <v/2H) 

r(N  +  2) 

(1  -7)N  (1  +  N7)  +  N(N  +  1) ^ (i  " |) 


N  - 1 


F(1-N-!;r-n7> 


1  <7  <2 


0<7  <  1 


f 


(34) 


10 


TR  4683 


where  y  is  the  scale  factor  in  EA  GLR  test  (24),  and  F  is  the  Gaussian  hyper- 
geometric  function.  The  factor  y  need  take  values  only  in  the  range  (0,  M); 
see  appendix  C  for  the  argument  for  M  =  2. 


Since,  for  general  M  and  N  ,  Pd  and  Pp  can  not  be  analytically  evalu¬ 
ated  for  the  EA  processor,  a  simulation  approach  to  determine  performance 
was  adopted.  Specifically,  in  order  to  determine  y  in  GLR  test  (24),  which 
realizes  specified  Pp  values  of  .  01  and  .  001,  40, 000  independent  trials  were 
conducted  for  the  EA  processor.  The  method  of  estimating  the  .  99  and  .  999 
quantiles  is  discussed  in  appendix  D. 


Since  the  determination  of  y  for  a  specified  Pp  was  accomplished  by 
a  simulation  based  on  40,  000  independent  trials,  the  estimate  of  y  has  a  non¬ 
zero  standard  deviation,  which  was  also  estimated  by  simulation.  In  the  next 
section,  the  detection  characteristics  for  the  EA  processor  are  associated  with 
the  estimated  y  ,  and  are  themselves  based  on  10,000  independent  trials. 
Therefore,  in  order  to  determine  the  reliability  of  the  plotted  curves,  PD  has 
also  been  determined  when  the  estimated  y  is  changed  by  + 1  standard  devi¬ 
ation,  as  determined  above.  Additionally,  the  standard  deviation  of  the  estimate 
of  Pq  ,  as  a  result  of  statistical  fluctuations  based  upon  10,  000  independent 
trials,  has  been  added.  Specifically,  when  the  true  value  of  the  detection 
probability  is  Pd  ,  the  standard  deviation  of  the  simulation  estimate  is 
(?d(1  “  Pd)/10*  000)1/' \  The  resultant  uncertainty  is  indicated  in  the  E A 
processor  performance  plots  by  broken  lines  parallel  to  the  central  curve.  We 
may  note  that  the  detection  characteristics  need  not  be  confined  to  the  uncer¬ 
tainty  region  described;  rather,  this  ribbon  describes  where  Pp  lies  with  fair 
confidence,  though  it  may  actually  deviate  beyond  the  ribbon. 


RESULTS 

Required  values  of  scale  factor  r  in  UA  test  (19)  are  presented  in 
table  B-l  for  Pp  =  10-rL  n  =  2,  3,4,  6,  8,  when  M  and  N  are  varied  as  follows: 
M  =  1,2, 3, 4,  5, 10  and  N  =  1,2, 3, 4,  6, 10,20,  «  .  The  performance  of  the  UA 
processor  is  shown  in  figures  1-30,  where  PD  is  plotted  versus  d-j-  (which 
is  a  voltage  measure  of  the  total  signal -to-noise  ratio),  with  Pp ,  M,  and 
N  as  parameters.  Each  figure  in  this  sequence  depicts  performance  over 
the  range  of  N,  indicated  above,  for  a  given  Pp  and  M.  The  entire  se¬ 
quence  of  figures  illustrates  performance  for  Pp  =  10-n,  n  =  2,  3,  4,  6,  8,  and 
M  =  1,2, 3,4,  5, 10.  The  curves  for  N  =  »  (known-noise  level)  have  also  been 
included  for  comparison. 
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Over  the  range  of  PF  considered,  UA  processor  performance  degrades 
sharply  when  N  is  too  few  (typically,  fewer  than  about  10  noise -alone  samples). 
Because  of  the  extremely  high  threshold  scale  factors  r  required  to  ensure 
low  PF ,  this  degradation  becomes  increasingly  severe  as  PF  is  lowered.  This 
is  illustrated  for  the  UA  processor  with  N  =  1 ,  where  the  required  threshold 
(scale  factor)  r  is  well  approximated  by 


1  M  +  1 

P„  "  2M 
F 


(N  =  l) 


(35) 


when  Pp  «  1 .  Thus  small  values  of  PF  require  exceedingly  large  r .  For 
example,  PF  =  10"^  requires  r  approximately  equal  to  10  .  (See  table  B-l 
for  other  values  of  N. ) 

The  improvement  attainable  in  UA  processor  performance  by  increasing 
N  is  more  obvious  in  figures  31-36.  These  figures  are  plots  of  required  dp 
in  decibels  (that  is,  10  log^d^,^)  )  to  realize  Pd  =  .  5,  with  Pp  and  M  as 
parameters,  versus  N.  The  horizontal  tic  marks  at  the  right  edge  of  these 
figures  show  the  required  value  of  dip  in  decibels  for  N  =  °°  and  are,  therefore, 
the  asymptotes  of  the  curves.  These  figures  also  indicate  that,  to  achieve  a 
given  PD  for  a  specified  dp  ,  N  must  be  increased  as  Pp  is  lowered. 

If  the  required  Pd  is  increased,  the  necessary  values  of  signal-to-noise 
ratio  dx  in  figures  31-36  would  increase  also.  However,  to  a  first  approxi¬ 
mation,  the  decibel  difference  between  the  values  of  dp  required  at  finite 
N  versus  those  required  at  N  =  *>  are  relatively  independent  of  Pd  ,  at  least 
for  N  >  6.  Thus,  for  a  range  of  Pd  including  approximately  (.5,  .99),  the 
additional  signal-to-noise  ratio  required  at  moderate  values  of  N  versus  those 
required  at  N  = 00  depends  only  on  N,  M,  and  Pp  .  If  the  curves  in  figures 
1-30  were  straight  lines  emanating  from  a  common  intersection  on  the  ordinate, 
this  behavior  would  be  precisely  correct.  Although  this  is  not  true  for  figures 
1-30,  an  extended  region  of  Pd  ,  where  extrapolations  of  the  curves  back  to 
dx  =  0  approximately  satisfy  this  requirement,  is  evident;  in  this  region,  the 
decibel  difference  discussed  above  is  relatively  independent  of  Pd  .  No  other 
independencies  of  N,  M,  or  P^,  have  been  observed. 

It  is  apparent  from  figures  1-36  that  UA  processor  performance 
steadily  degrades  as  the  signal  energy  fractionalization  increases,  that  is,  as 
M  increases,  for  fixed  dx .  For  example,  for  PF  =  10  "2«  PD  =  .  5,  and 
N  =  2,  the  required  value  of  dx  is  10.  4  for  M  =  5  versus  5.4  for  M  =  1, 
a  5.  7-dB  increase.  The  situation  is  somewhat  alleviated  for  larger  N  ;  for 
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-2 

example,  at  Pp  =  10  ,  Pp)  =  .  5,  and  N  =  10,  the  required  value  of  d-p  is 

4.  8  for  M  =  5  versus  3.  2  for  M  =  1,  only  a  3.  5-dB  increase.  Also  the  degra¬ 
dation  is  not  quite  so  severe  for  smaller  Pp ;  for  example,  at  Pp  =  10  , 

PD  =  .  5,  and  N  =  20,  the  required  values  of  dT  for  M  =  5  versus  M  =  1  differ 
by  3. 1  dB. 

When  comparing  the  results  of  this  report  with  those  in  reference  1,  note 
that  a  direct  overlay  does  not  constitute  a  valid  comparison.  For  example,* 
figure  17  in  reference  1,  when  overlaid  upon  figures  31-36  here,  yields  curves 
that  intersect  for  some  values  of  N .  The  reason  is  that  N  ,  when  specified 
here,  corresponds  to  2N  degrees  of  freedom  in  the  random  variable  on  the 
right  side  of  test  (19);  on  the  other  hand,  the  variable  N  in  reference  1  corre¬ 
sponds  to  only  N  degrees  of  freedom.  Before  a  valid  comparison  is  drawn, 
this  factor  of  2  difference  in  the  interpretation  of  N  in  the  two  reports  must  be 
accounted  for.  This  observation  is  also  consistent  with  the  fact  that  the 
envelope  of  the  output  of  a  narrowband  filter  need  be  sampled  only  half  as  often 
as  a  low-pass  process  with  the  same  (positive-frequency)  bandwidth  as  that  of 
the  narrowband  filter.  If  the  potential-signal  (matched)  filter  outputs  are  each 
sampled  once,  as  indicated  in  (5),  only  this  correction  is  needed  in  translating 
results.  It  will  then  be  found  that  the  curves  in  reference  1  lie  everywhere 
below  those  of  figures  31-36;  that  is,  phase-coherent  operation  yields  better 
performance  than  phase-incoherent  operation,  as  is  to  be  expected. 

Estimated  values  of  the  threshold  y  for  the  EA  processor  are  given  in 
table  D-l.  As  noted  in  the  previous  section,  the  values  must  be  estimated  by 
simulation  and,  thus,  are  not  exact  (except  for  M  =  2,  which  is  treated  in 
appendix  C). 

Simulated  performance  of  the  EA  processor  is  shown  in  figures  37-46; 
there  Pd  has  been  plotted  versus  dx  for  Pp  =  10“2  and  10“3,  with  M  and  N 
as  parameters.  It  should  be  emphasized  that,  for  a  given  dp  ,  all  the  M  signal 
amplitudes  were  taken  to  be  equal  in  the  simulation.  The  performance  curves 
for  M  =  1  have  not  been  presented  because,  for  this  case,  the  UA  and  EA 
processors  are  equivalent;  for  M  =  1,  figures  1-5  define  the  performance  of 
the  UA  and  EA  processors.  Figures  37-46  present  EA  processor  perform¬ 
ance  for  M  =  2,3,4,5,10  and  N  =  1, 2, 3,4,  6, 10,  20.  (An  analytical  method 
of  determining  performance  for  N  =  »  is  described  in  appendix  D;  it  was  not 
pursued  here,  however. ) 


♦Notice  that  the  decibel  definition  of  signal-to-noise  ratio  in  reference  1, 
eq.  (17)  et  seq.  ,  differs  from  that  adopted  here. 
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Consistently,  EA  processor  performance  is  superior  to  that  of  the  UA 
processor.  *  It  is  to  be  noted  that  the  smaller  the  N,  the  greater  the  superiority 
of  the  EA  processor.  The  principal  reason  for  this  performance  difference 
(other  than  the  equal  signal  amplitudes  used  in  the  simulation)  can  be  seen 
in  the  forms  of  the  EA  and  UA  tests,  (24)  and  (19),  respectively.  The  EA  test 
(except  for  M  =  1)  also  uses  the  potential-signal  samples  being  tested  to  help 
form  a  better  noise-level  estimate,  that  is,  the  right  side  of  (24);  the  UA  proc¬ 
essor  restricts  itself  by  using  noise-alone  samples  to  form  a  noise-level 
estimate,  that  is,  the  right  side  of  (19).  For  small  N,  the  use  of  M  +  N  samples 
significantly  improves  the  noise-level  estimates.  As  N  is  made  larger, 
the  potential-signal  samples  become  less  important  in  estimating  the  noise 
level. 


Another  example  of  this  effect  is  evident  from  the  performance  of  the  EA 
processor  for  fixed  Pp  and  dq> ,  say  10-2  and  10,  respectively.  Then,  for 
N  =  0,  as  M  increases  from  2  to  4,  for  example,  Pj)  improves  from  .  04  to 
.  18.  For  N  =  1,  the  change  in  PD  is  from  .53  to .  62;  for  N  =  2,  Pjj  decreases 
from  .  95  to  .  925;  and  for  N  >  2,  performance  degrades  as  M  increases  from 
2  to  4.  The  reason  for  this  effect  is  that  for  large  N,  where  a  good  estimate 
of  noise  power  results,  fractionalization  of  a  fixed  total  amount  of  signal  energy 
into  more  channels  hampers  EA  processor  performance  because  the  recombina¬ 
tion  of  signal  energy  is  incoherent.  However,  for  N  =  0,  where  a  noise-alone 
power  estimate  is  not  available,  fractionalization  of  a  fixed  total-signal  energy 
into  more  channels  affords  some  (albeit  poor)  estimate  of  the  noise  power, 
thereby  establishing  a  threshold  for  specified  Pp .  (Probably  if  too  large  an  M 
is  chosen,  the  fractionalization  will  be  too  great  and  performance  will  degrade 
with  increasing  M. )  For  intermediate  N,  such  as  N  =  1,  a  transition  be¬ 
tween  these  two  effects  takes  place. 

The  improvement  in  performance  of  the  EA  processor,  as  the  number  of 
noise-alone  samples  is  increased,  is  shown  explicitly  in  figures  47-51,  where 
the  dj,  (in  decibels)  required  to  realize  Pjj  =  .  5  versus  N,  with  Pp  and  M 
as  parameters,  is  plotted.  The  curve  for  M  =  1  is  not  presented  since  it  is 
identical  to  figure  31,  as  noted  earlier.  The  comments  given  above  concerning 
the  use  of  potential-signal  samples  on  both  sides  of  the  GLR  test  are  illustrated 
in  figures  50  and  51.  In  these  plots,  for  N  significantly  smaller  than  M,  the 
leveling  off  of  the  curves  shows  that  the  potential-signal  samples  in  these 
instances  are  more  important  than  the  noise-alone  samples  in  establishing 
stable  noise-level  estimates. 


*We  reiterate  that  the  signal  amplitudes  have  been  set  equal  to  each  other 
in  the  simulation  for  the  EA  processor.  This  will  be  discussed  again  later. 
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DISCUSSION 

Operating  characteristics  for  detecting  phase-incoherent  signals  in 
Gaussian  noise  of  unknown  level  have  been  presented.  The  two  processors  con¬ 
sidered  were  ML  detectors  realized  under  different  assumptions  about  signal 
amplitudes.  In  the  first  case,  the  signal  amplitudes  were  assumed  unequal  and 
unknown  from  sample  to  sample.  In  the  second  case,  the  signal  amplitudes 
were  assumed  equal,  though  unknown.  Since  the  GLR  tests  contain  no  claims 
to  optimality,  it  may  be  questioned  whether  other  processors  obtained  from 
different  and,  presumably,  ad  hoc  tests  might  outperform  the  processors  con¬ 
sidered  here.  In  order  to  partially  answer  this  question,  a  number  of  other 
tests,  in  addition  to  the  UA  and  EA  GLR  tests,  were  compared  by  means  of 
simulation.  We  shall  now  very  briefly  consider  some  reasonable  ad  hoc  tests 
and  the  results  of  the  simulations  of  their  performance. 

To  facilitate  the  discussion,  we  define  the  following  quantities  in  terms 
of  previously  used  notation: 


(36) 


Recall  that  (z^J  and  |zj  |  are  envelope  samples  at  the  outputs  of  the  potential- 
signal  and  noise-alone  filters,  respectively.  Using  the  above  definitions,  we 
can  write  the  total-sample  variance  as 


Var 


2  ~ 
S -  Mm  +  S 


Nm 


M  +  N  -  2 


(37) 


Then  the  UA  test  in  the  new  notation  can  be  written  as 


S  >  S 
M  <  C  N 


(38) 
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where  c  (>  0)  is  a  scaling  constant,  and  the  EA  test  can  be  written  as 


(39) 


The  additional  ad  hoc  tests  considered  here  are  indicated  in  table  1.  All  the 
tests  considered  are  invariant  to  level  scaling;  that  is,  they  yield  CFAR 
processors. 


Table  1.  Examples  of  Ad  Hoc  Processors  Invariant  to  Level  Scaling 


Case 

Processing 

(a) 

m  <  c  (Var)1^2 

(b) 

~  *  __  1/2 

m  -  m  £  c  (V ar) 

(c) 

~  *  /s  +s\1/2 

m  -  m  ^  c^m  +  n  j 

<d) 

m  >  c  (S/N)1/2 

(e) 

m  ^  cm 

CO 

m  -  m  ^  c  (S/N)1^2 

(g) 

Mm  - Nm  5  C(M  +  n) 

<h) 

Mm  -  Nm  £  c  (S/N)1^2 

Note:  If  the  upper  inequality  is  satisfied,  is  declared  true. 

If  the  lower  inequality  is  satisfied,  HQ  is  declared  true. 
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Case  (b)  in  table  1  is  an  estimate  of  the  usual  deflection  criterion  (refer¬ 
ence  5)  compared  with  threshold  c.  Cases  (c)  and  (d)  are  modifications  of  the 
EA  processor,  and,  analogously,  case  (a)  may  be  thought  of  as  a  modification 
of  the  deflection  criterion.  Case  (f)  is  suggested  by  a  comparison  of  cases  (c) 
and  (d).  Cases  (g)  and  (h)  are  extensions  of  cases  (c)  and  (f),  where  the  sample 
means  have  been  weighted  according  to  the  number  of  samples  used.  Finally, 
case  (e)  is  a  straightforward  comparison  of  sample  means. 


The  processors  in  table  1  and  the  EA  and  UA  processors  were  compared 
by  simulating  their  performance  over  a  wide  range  of  values  of  M,  N,  and  dx . 
When  equal  signal  amplitudes  were  used,  the  EA  processor  outperformed  all 
others.  For  unequal  received  signal  amplitudes,  it  has  already  been  noted  that 
the  performance  of  the  UA  processor  is  unaffected  as  long  as  dx  is  kept  con¬ 
stant.  In  fact,  the  UA  processor  performed  best  for  some  examples  of  unequal 
signal  energy  fractionalization  used  in  the  simulation,  such  as  0  in  one  of  the 
M  potential-signal  branches.  As  an  example,  the  EA  processor  for  M  =  4  was 
subjected  to  unequal  signal  amplitudes,  namely,  zero  strength  in  one  channel 
and  equal  signal  strength  in  the  remaining  channels.  (A  related  situation  was 
treated  in  reference  2. )  The  detection  characteristics  for  Pp  =  10-2  and  10-3 
are  given  in  figures  52  and  53.  Comparison  with  figures  41  and  42  shows  the 
degradation  to  be  significant,  especially  for  small  values  of  N .  Thus  the  EA 
processor  is  sensitive  to  the  exact  received  signal  energy  fractionalization 
amongst  branches,  whereas  the  UA  processor  is  completely  insensitive  to  this 
effect.  In  all  the  cases  studied,  one  of  the  two  ML  processors  (EA  or  UA)  out¬ 
performed  all  the  other  processors  in  table  1. 


In  the  foregoing  discussion,  we  considered  the  EA  processor  per¬ 
formance  when  unequal  signal  energy  fractionalization  occurred  in  the  M 
potential-signal  branches.  We  have  also  noted  that  UA  processor  performance 
is  unaffected  by  the  manner  in  which  the  received  signal  energy  is  fractionalized. 
Alternatively,  if  {d^}  are  considered  to  be  random  variables  with  known  PDFs, 
the  present  results  can  be  averaged  over  the  given  PDFs  to  obtain  the  average 
P£>  for  fluctuating  signal  amplitudes.  For  the  EA  processor,  graphical  inte¬ 
gration  is  required  for  this  averaging  procedure;  for  the  UA  processor, 
analytical  evaluation  is  possible  for  some  PDFs.  The  case  of  dx/2  having  a 
chi-squared  PDF  will  be  considered  in  part  m  of  this  sequence  of  technical 
reports  (reference  1). 
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When  the  threshold  y  is  varied  by  +1  standard  deviation  away  from  its 
estimated  value  for  a  specified  Pj? ,  large  percentage  changes  in  Pj)  for  the 
EA  processor  sometimes  occur.  This  is  particularly  true  for  small  N  and  Pp . 
In  these  cases,  the  EA  processor  is  very  sensitive  to  its  precise  threshold 
setting  and  is  not  a  viable  processor.  The  UA  processor  is  less  sensitive  to 
threshold  setting  because  the  M  potential-signal  samples  are  not  in-bred;  that 
is,  they  are  used  only  on  the  left  side  of  the  GLF  test.  However,  no  analysis 
of  the  sensitivity  to  the  threshold  was  conducted  for  either  processor. 
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Figure  1.  Detection  Characteristics  for  UA  Processor;  M  =  1,  Pp  =  10 
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Figure  2.  Detection  Characteristics  for  UA  Processor;  M  =  1,  PF  =  10  3 
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Figure  3.  Detection  Characteristics  for  UA  Processor;  M  =  1,  Pp  =  10 ~4 
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Figure  4.  Detection  Characteristics  for  UA  Processor;  M  =  1,  PF  =  10-6 
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Figure  5.  Detection  Characteristics  for  UA  Processor;  M  =  1,  Pp  =  10-P 
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Figure  6.  Detection  Characteristics  for  UA  Processor;  M  =  2,  Pp  =  10“2 
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Figure  7.  Detection  Characteristics  for  UA  Processor; 
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Figure  8.  Detection  Characteristics  for  UA  Processor;  M  =  2,  Pp  =  10“4 
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Figure  9.  Detection  Characteristics  for  UA  Processor;  M  =  2,  Pp  =  10-6 
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Figure  10.  Detection  Characteristics  for  UA  Processor;  M  =  2,  Pp  =  10  8 
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Figure  11.  Detection  Characteristics  for  UA  Processor;  M  =  3,  Pp  =  10 
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Figure  12.  Detection  Characteristics  for  UA  Processor;  M  =  3,  Pp  =  10~3 
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Figure  13.  Detection  Characteristics  for  UA  Processor;  M  =  3,  PF  =  10"^ 
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Figure  14.  Detection  Characteristics  for  UA  Processor;  M  =  3,  PF  =  10  6 
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Figure  15.  Detection  Characteristics  for  UA  Processor;  M  =  3,  PF  =  10-8 
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Figure  16.  Detection  Characteristics  for  UA  Processor;  M  =  4,  Pj,  =  10  2 


34 


TR  4683 


dT 

Figure  17.  Detection  Characteristics  for  UA  Processor;  M  =  4,  Pp  =  10"^ 
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Figure  18.  Detection  Characteristics  for  UA  Processor;  M  =  4,  Pp,  =  10  ^ 
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Figure  19.  Detection  Characteristics  for  UA  Processor;  M  =  4,  Pp  =  10“® 
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Figure  20.  Detection  Characteristics  for  UA  Processor;  M  =  4,  PF  =  10"® 
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Figure  21.  Detection  Characteristics  for  UA  Processor;  M  =  5,  PF  =  10 
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Figure  22.  Detection  Characteristics  for  UA  Processor;  M  =  5,  Pj,  =  10“^ 
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Figure  23.  Detection  Characteristics  for  UA  Processor;  M  =  5,  PF  =  10  4 
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Figure  24.  Detection  Characteristics  for  UA  Processor;  M  =  5,  Pj,  =  10-6 
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Figure  25.  Detection  Characteristics  for  UA  Processor;  M  =  5,  P-p  =  10-8 
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Figure  26.  Detection  Characteristics  for  UA  Processor;  M  =  10,  Pj,  =  10-2 
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Figure  27.  Detection  Characteristics  for  UA  Processor;  M  =  10,  Pp  =  10"^ 


45 


TR  4683 


46 


TR  4683 


dT 


Figure  29.  Detection  Characteristics  for  UA  Processor;  M  =  10,  Pp  =  10  ® 
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Figure  30.  Detection  Characteristics  for  UA  Processor;  M  =  10,  P„  =  10"® 
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Figure  31.  Required  Signal-to-Noise  Ratio  for  UA  Processor; 
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Figure  32.  Required  Signal-to-Noise  Ratio  for  UA  Processor; 
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Figure  33.  Required  Signal-to-Noise  Ratio  for  UA  Processor; 
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Figure  34.  Required  Signal -to -Noise  Ratio  for  UA  Processor; 
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Figure  35.  Required  Signal-to-Noise  Ratio  for  UA  Processor; 
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Figure  36.  Required  Signal-to-Noise  Ratio  for  UA  Processor; 
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Figure  37.  Detection  Characteristics  for  EA  Processor;  M  =  2,  Pp  =  10-^ 
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Figure  38.  Detection  Characteristics  for  EA  Processor;  M  =  2,  PF 
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Figure  39.  Detection  Characteristics  for  EA  Processor;  M  =  3,  Pp  =  10 
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Figure  40.  Detection  Characteristics  for  EA  Processor;  M  =  3,  Pp 
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Figure  41.  Detection  Characteristics  for  EA  Processor;  M  =  4,  Pp  =  10  2 
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Figure  42.  Detection  Characteristics  for  EA  Processor;  M  =  4,  Pp  =  10“3 
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Figure  43.  Detection  Characteristics  for  EA  Processor;  M  =  5,  Pp  =  10  ^ 


61 


TR  4683 


Figure  44. 


Detection  Characteristics  for  EA  Processor;  M  =  5,  Pp  =  10“^ 
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Figure  46.  Detection  Characteristics  for  EA  Processor;  M  =  10,  Pp  =  10  ^ 
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Figure  47.  Required  Signal-to-Noise  Ratio  for  EA  Processor;  M  =  2 
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Figure  48.  Required  Signal-to-Noise  Ratio  for  EA  Processor;  M  =  3 
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Figure  49.  Required  Signal-to-Noise  Ratio  for  EA  Processor;  M  =  4 
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Figure  50.  Required  Signal-to-Noise  Ratio  for  EA  Processor;  M  =  5 
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Figure  51.  Required  Signal-to-Noise  Ratio  for  EA  Processor;  M  =  10 
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Figure  52.  Detection  Characteristics  for  EA  Processor 
with  Unequal  Amplitudes;  M  =  4,  P^,  =  10 
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Figure  53.  Detection  Characteristics  for  EA  Processor 
with  Unequal  Amplitudes;  M  =  4,  PF  =  10~3 
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Appendix  A 

PREDETECTION  PROCESSING  AND  STATISTICS 


Initial  receiver  processing  consists  of  computing  the  complex  quantities 

J  dt  Fk(t)  Vk(t)  =  zk  ,  1  <  k  <  M  ,  -  (A-l) 

for  the  potential-signal-plus-noise  branches,  and 

J  dt  F  (t)  V  (t)  =  z  ,  1  <  j  <  N  ,  (A-2) 

for  the  noise-alone  branches.  (Integration  is  over  the  range  of  nonzero  inte¬ 
grand. )  In  the  above  expressions,  (V^t)}  and  (Vj(t)}  are  the  complex 
envelopes  of  the  received  processes,  and  (F^t)}  and  (Fj(t)}  are  receiver 
filter  complex  envelopes  that  may  or  may  not  be  optimum.  *  The  complex 
random  variables  and  {zj}  are. Gaussian  because  the  input  noises  are 

stationary  Gaussian  processes. 


Denoting  ensemble  averages  on  the  noise  by  overbars,  we  can  write 

N  (t)  =  0  ,  N.(t)  =  0  ,  all  k ,  j  .  (A-3) 

K  J 

In  (A-3),  {Nj^t)}  and  (Nj(t)}  are  the  complex  envelopes  of  the  received  noise 
waveforms.  Then  the  noise  average 


forHl 

0 ,  for  H 

o 


1  <  k  <  M  , 


(A -4) 


where  {Sk(t)}  are  the  complex  envelopes  of  the  actual  received  signal  wave¬ 
forms.  Similarly 

z  =  0  for  H  and  H  ,  1  <  j  <  N  .  (A-5) 

j  1  o  -  - 


*For  example,  mismatch  or  misalignment  in  frequency  and/or  time  are 
allowed. 
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Noting  that 

zk  "  \  =  \  dt  Ffc(t)  Nk(t)  ,  1  <  k  <  M  ,  (A-6) 

for  H-^  and  HQ ,  it  is  seen  that 


(zk  ‘  zk)(zj!  -  *j)  =  0  for  aU  k' 1  ’  <A‘7> 

since,  for  complex  envelopes  of  stationary  processes, 


Nk(t)  (u)  =  0  for  all  k,  J2  .  (A-8) 

Also,  writing  an  element  of  the  received  noise  crosscorrelation  matrix  as 


Nk(t)  N*(u)  =Rki(t-u)  ,  (A-9) 

we  obtain 

(zk_ik)(zi! •z«)*=ndtd"Fk(*)Fi(u)Rki(t_u)52<,ki-  <a-io) 

Similarly,  for  the  noise-alone  branches, 


z.  z.  =0  (A- 11) 

j  i 

and 


z.  zj  =  dt  du  Fj(t)  F*  (u)  S  (t  -  u)  =  27Tj  .  (A- 12) 


In  the  present  investigation,  attention  is  restricted  to  the  situation 


It  is  also  assumed  that 


2  . 

-  'o  SU  ’ 


2 

=  <r  S.A  . 
o  JJ2. 


(A-13) 


^zk  -  zk^  z*  =  0,  all  k,  a  ,  (A- 14) 

which  corresponds  to  uncorrelated  noises  in  the  potential-signal  branches  and 
noise-alone  branches.  An  example  of  (A-13)  and  (A-14)  is  afforded  by  disjoint 
equal -gain  filters  and  a  fairly  flat  noise  spectrum  over  the  complete  receiver 
filter  bank. 
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Appendix  B 

PERFORMANCE  EVALUATION  OF  UA  PROCESSOR 


Let 


z  =  B  exp(i0 ),  1  <  k  <  M;  =  B.  exp(i<T.j,  1  <  j  <  N  ,  (B-l) 

k  k  k  J  J  J 


where  {z^}  and  {  z'j}  are  the  filter  outputs  defined  previously  (cf.  appendix  A). 
If  we  let  vector 

B  =  ^  .  Bm  B^  .  . .  BnJ  ,  (B-2) 


then,  if  the  Gaussian  character  of  {z^}  and  (“zj}  and  the  statistics  presented  in 
appendix  A  are  used,  the  PDF  of  B  under  hypothesis  is 


P(B)  = 


M 

XT 

r  \ 

’  Bk  ‘  2Ako  Bk  cosWko  "  V  +  Nm] 

11 

J  dek  2  exp 

2 

k=l 

J2*  27t<to 

1 - 

to 

o 

1 _ 

N  b. 

n  \  a5s  “Texp 

’=1  *2ir  2na 


1=1 

M 

n 

k=l 


B, 


B2"1 
J 


2<x 

o  J 


(B-3) 


exp 


2  2  n 

Bk  +  Ako 


2  a 


I  f^”8' 


N 

n 

j=i 


B. 


exp 


B21 
J 

2.2 

0J 


B,  >  0,  B.  >  0  . 
k  J 


Here  A^q  is  the  actual  amplitude  of  the  complex  sample  of  the  signal  output  of 
the  filter  F^t)  ;  that  is, 

(B-4) 


ko 


=  |J  dtFk(t)Sk(t) 


1  <  k  <  M  , 


where  Sk(t)  is  the  actual  received-signal  complex  envelope.  The  parameter  a £ 
in  the  PDF  (B-3)  is  proportional  to  the  actual  output  noise  variance  of  a  typical 
filter  (see  (A-10)  and  (A-13)): 


?o=\\\  dt  du  Fx(t)  F*(u)  Rn(t  -  u)  . 


(B-5) 
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The  GLR  test  of  interest  is  obtained  from  (19)  and  (B-l): 

M  N 


j=l 


(B-6) 


k=l 


Let 


M  /B. 


U 


sE 


XT  /C?  v2 

N  /B 


k=i  '  ° 

Then  test  (B-6)  can  be  written  as 


j=l  °' 


''Sco 


dk=  . 


(B-7) 


U  %  «V  . 


(B-8) 


The  random  variables  U  and  V  are  statistically  independent  of  each  other;  see 
(B-3).  Therefore  the  characteristic  function  of  U  is  (reference  6,  6.  631  4) 


ty*)  =  exp(i{U)  = 


M 

■n|( 


dx  exp(i£x  )  x  exp 


[- 


2  ,2n 

x  +dk 


I  (d,  x) 
o'  k  1 


-M  /‘dT( 
=  (l-i2()  M  exp 


where 


M 


1  -i2  £/  ’ 


M 


K  5E  K=Et 


i/dtFk<t)sk(t>r/2 


K  •  fe  j//d‘  du  Fjffl  F»(u)  Rn(t-u) 


From  (B-9),  the  PDF  of  U  can  be  determined  to  be 

2  . M_1 

x  +  d 


(B-9) 


(B-10) 


(B-ll) 


In  a  similar  manner,  the  PDF  of  V  can  be  shown  to  be 

P  (y)^  exp  ('-X)  iY/2>N~L  ,  y>0. 

Hyvj'/  2  vy  2/  »  y  *  ” 


(B-12) 
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From(B-8),  (B-ll),  and(B-12),  there  is  obtained  for  the  probability  of 
detection 


PD  =  Prob  (U  >  aV) 


dx  pu(x) 


(y)  • 


(B-13) 


Consider  the  inner  integral  in  (B-13): 


(y) 


(y/2)N-1 

(N  -  1) ! 


=  1  -  exp 


N-l 


(-rJE  si 


k=0 


(B-14) 


Substitution  of  (B-ll)  and  (B-14)  in  (B-13)  and  a  change  of  variable,  x  =  t2  ,  yields 

r..“  -[■  fr-aiwv- 

d„  wa  •'0 


P„  =  1 

‘S?  k=0 

Recall  the  identity 


*M-l^dT  ^  .M-l  JM-l^ldT^ 


(B-16) 


then  the  integral  in  (B-15)  can  be  evaluated  from  reference  6,  6.  631  1  as 


r  tM+2k  ft2/  nl 
1 dt  “iTT  expL'  Tl1  -  «)J  JM-l(IdTt) 


M-l 

r(M+k) 

M/a+l\M  +  k 

2  br)  (n-1)! 


(B-17) 


=  oV.  a  \M+k  jM_1  (M  +  k  -  1) !  /jr  _a\  /  a  \ 

2  \a  +  l)  dT  (M-l)!  6X15 \ 2  a  +  l/lFl\  k,M,“  2  a+1/’ 
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where  1F^(a;bpt)  denotes  the  confluent  hypergeometric  function.  One  form  for 
the  probability  of  detection  for  the  UA  processor  is  thus 


The  probability  of  false  alarm  is  obtained  when  d>p  =  0  in  (B-18): 


(B-20) 


An  alternative  form  for  Pjj  can  be  obtained  by  first  noting  that  (refer¬ 
ence  7,  13.6.9) 


'M  -  1  +  k\ 


1^'1  (-k;  M;x)  =L^M-1)(x) 


(B-21) 


where  the  generalized  Laguerre  polynomials  (l^  ^(x))  satisfy  the 


recurrence 


relations 


(M-l) 
Jo 


L'*  *  (x)  =  1,  L'~  ~  (x)  =  M  -  x  , 


(M-l) 
J1 


Substitution  of  (B-21)  in  (B-18)  gives 

/  dt  \  M  N-l 

pD=1‘exp(”r  E-—kLk 

'  k=0<a+1> 


(M-l)  (  dT  a 
2  a  +  1 


(  iT  \  M  N-l  /  d  \ 

S1-“p(- T  STl)(sTl)  ZX(‘T  M)- 

k=0 


(B-22) 


if-1* (X,  =  2k-2;M-* iff (X) - if;1* k > 2 . 


(B-23) 
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(For  M  =  1,  the  upper  line  of  (B-23)  agrees  with  reference  2,  eq.  (14).)  It  is 
seen  that  (B-23)  is  a  series  of  positive  terms  that  satisfy  the  recursion  relations: 

]Uf  —  V 

G(x,a,M)  =  l,  G  (X,ar,  M)  = 


a  +  1  * 


GJx,a,  M)  =  ~T7~t  A~X  Gwn  (x,a,  M)  -  — —  ^  G^x.a,  M),  k>  2  . 


(B-24) 


k(ar  + 1)  k-1 


k(a+l) 

We  consider  now  an  alternative  expression  for  the  probability  of  detection 
for  the  UA  processor; 


PD  =  Prob  (U  >  crV) 


dy  Pv(y) 


P„W  . 


(B-25) 


Using  (B-ll),  we  can  write  the  inner  integral  in  (B-25)  as 


where  Qj^(a,  b)  is  the  generalization  of  the  Q-function  (reference3).  Then 
Pq  ,  (B-25),  becomes,  with  the  help  of  (B-12), 
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which  by  a  change  of  variable,  y  =  t2,  can  be  written  in  the  form 

i.2\  t2(N-l) 

2N_1(N  -  1)! 


r  t\  t';  r 

P  =\  dttexp^-— j— - |Q(d 

U  •'0  x  z/  2w  fN  -  1\!  L  1 


•  i/2  M 
a  t) 


at2  *  4\  W  Ll/Z  \k 

+ “p' — — )  22  hr  v  °  •>] 


ri v  ‘‘t 


(B-28) 


EI1  +  I2- 


Now  consider  the  term  1^  ,  which  is  defined  as 


■x-I 


2  t2N_1 

I, =  \  dt  exp(-t  /2)  nrr - 

J0  2N_1(N-1)! 


Q(dT,  «1/2t) 


-r 


(2  2  \ 

X  +0 

- o - )  I  (dmX)  • 


(B-29) 


2  /  o'  T' 


By  interchanging  the  order  of  integration  in  (B-29)  and  then  making  a  change  of 
variable,  v  =  t2/2,  we  obtain 


’■■1 


I  =  |  dx  x  expl- 


2  ^2\ 
x  +<y 


x/al/2 

f  2 

IQ(dT  x)  J0  dt  exP(“t  /2) 


2N-1 


2  (N  -  1) ! 


=  1  - 


&£ 

k=0  ** 


2k+l 


dx  r-  exp  (-  — 5^)  ID(dT  X)  exp  (-  fj)  , 


(B-30) 


(2«) 
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which  can  be  evaluated  as  (reference  6,  6.  631  1) 


d2\N-l 


1=1  — expl — — )\  ^  — - — -  F 

“  V*  yK(«+D  1 


k+1-'l'fa-7T 


=  1  "^TT  6XP 


(  d2  X  -N-1  /  d2 

v  2  a+l/s<“+i)k  1  A  ’  ’  2  a+ 


(B-31) 


N-l 


with 


=  1  ”  a+l  “p(-  T  Tk  Lk(-  T  srij ' 

k=0'a+i) 


L(0)  (x)  =  L  (x) 
n  '  n 


(B-32) 


The  remaining  term  1 2  in  (B-28)  can  be  developed  as  (reference  6,  6.  631  1) 

.2\ 


I2=exp 


/  d“  \  .  M-l  .  .  k/2  .00 

ldtt 


2N-l+k 


exp 


rT<1+“>)yV1/2‘> 


-^-(4  £T‘;“)fe)‘  ...(.-....‘.-i-u 


(1+a) 


k=l 


2  1  +  or  /• 

(B-33) 


Substituting  (B-31)  and  (B-33)  in  (B-28)  yields 


a 


PD  =  1  -  1+or 


/  d2  \  N-l  /  d2 

k=0  (1+ot) 


(B-34) 
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A  table  of  required  threshold  values  r  in  (B-6)  is  presented  in  table  B-l. 
This  table  was  determined  from  (B-35)  and  (B-38). 


Table  B-l.  Threshold  Values  r  for  UA  Processor 


O 

1 

to 

CO 

1 

o 

rH 

10  "4 

10-6 

00 

O 

iH 

1 

99. 

999. 

9999. 

999999. 

99999999. 

2 

18. 

61. 2455532 

198. 

1998. 

19998. 

3 

10. 9247665 

27. 

61. 6330407 

297. 

1389.47665 

4 

8. 64911064 

18.4936530 

36. 

122.491106 

396. 

6 

6.92660814 

12. 9736660 

21. 8495330 

54. 

123. 266081 

10 

5. 84893192 

9.95262315 

15. 1188643 

29.8107171 

53.  0957344 

20 

5. 17850824 

8.25075089 

11.  6978638 

19.  9052463 

30.2377286 

OO 

4. 60517019 

6. 90775528 

9. 21034037 

13. 8155106 

18.4206807 

h-1 

O 

1 

to 

op 

o 

rH 

10-4 

i<f6 

O 

00 

1 

99.2493719 

999. 249937 

9999.24999 

999999. 250 

99999999.  25 

2 

15. 9770249 

53.4358291 

171.870780 

1730. 71738 

17319. 1747 

3 

9. 14830103 

21.  9235414 

49.4186582 

236.  233970 

1103.33419 

4 

7. 00607662 

14. 3915845 

27.4928357 

92. 1698584 

296. 668412 

6 

5.41195143 

9.  63272610 

15. 7925977 

38.  0476483 

85.9572018 

10 

4.43069016 

7. 09603407 

10.4147233 

19.7920962 

34. 6109950 

20 

3. 82829355 

5. 69813414 

7. 75922632 

12.  6065444 

18.  6643142 

OO 

3.31917603 

4. 61670674 

5. 87818562 

8.  34421042 

10. 7678926 

X 

10-2 

io"3 

1 

o 

1—1 

io'6 

00 

o 

i-H 

1 

99.3325889 

999. 333259 

9999.33333 

999999.333 

99999999.3 

2 

15. 2068649 

50. 5250219 

162. 187065 

1631. 88194 

16328. 8205 

3 

8.46612534 

20.  0296547 

44. 9092656 

213. 942069 

998.499699 

4 

6.37068073 

12. 8580261 

24.3565038 

81. 1059802 

260. 530367 

6 

4. 82057350 

8.37881423 

13.5601910 

32.2606939 

72. 5053353 

10 

3. 87142682 

6. 01860847 

8.67889223 

16. 1732390 

27.9997019 

20 

3.29101239 

4. 73056833 

6.30305962 

9.97677756 

14. 5488339 

00 

2. 80198231 

3.74295742 

4. 64272354 

6.37638938 

8. 06043792 
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Table  B-l  (Cont'd).  Threshold  Values  r  for  UA  Processor 


\F 

h-1 

o 

1 

to 

H-1 

O 

1 

CO 

-4 

10 

CO 

1 

o 

T— 1 

00 

1 

o 

iH 

1 

99.3742148 

999.374922 

9999.37499 

999999.375 

99999999.4 

2 

14. 7988888 

48. 9961888 

157. 112691 

1580. 1387.1 

15810.3883 

3 

8. 10165137 

19. 0303331 

42.  5408682 

202. 267363 

943. 627975 

4 

6.02887011 

12. 0455412 

22. 7056389 

75.3107695 

241.  628275 

6 

4.49936528 

7.71035231 

12.3807118 

29.2278151 

65.4777610 

10 

3.56441205 

5.43999319 

7.75725282 

14.2739056 

24. 5490000 

20 

2. 99298087 

4.20703658 

5. 52574668 

8. 59348159 

12.4009667 

CO 

2.51127938 

3.26556019 

3. 97845350 

5. 33761427 

6.  64618478 

\F 

10-2 

CO 

1 

o 

iH 

10-4 

CD 

1 

o 

rH 

h-1 

o 

1 

00 

1 

99.3991960 

999.399920 

9999.39999 

999999.400 

99999999.4 

2 

14. 5459008 

48. 0525891 

153.984789 

1548. 25988 

15491. 0000 

3 

7. 87411853 

18. 4109248 

41.0768153 

195. 062504 

909.776375 

4 

5. 81429386 

11.5400561 

21.6825553 

71.7299629 

229. 958941 

6 

4. 29605440 

7.29202903 

11. 6467066 

27.3500642 

61. 1356414 

10 

3.36818639 

5. 07524621 

7. 18053924 

13. 0944804 

22. 4142485 

20 

2. 80054511 

3. 87438608 

5.03629308 

7.73123567 

11. 0695697 

00 

| 2.32092512 

2. 95882985 

3. 55640140 

4. 68630468 

5.76639616 

X 

<N 

1 

O 

iH 

I-* 

o 

1 

CO 

h-1 

O 

1 

CD 

1 

O 

iH 

00 

1 

o 

iH 

1 

99.4491708 

999.449917 

9999.44999 

999999.450 

199999999.45 

2 

14.0196087 

46. 1002576 

147.522750 

1482.43958 

14831.  5970 

3 

7.39583189 

17. 1201102 

38.  0358331 

180. 127861 

839. 636486 

4 

5. 35909494 

10.4796828 

19.5469671 

64. 2836270 

205.  719307 

6 

3. 85843310 

6.40480557 

10. 1014505 

23.4241458 

52.  0825303 

10 

2. 93773528 

4. 28996645 

5. 95161262 

10. 6085787 

17. 9392886 

20 

2.36887612 

3. 14498953 

3.97720089 

5. 89378598 

8. 25675959 

00 

1.  87831174 

2.26573733 

2. 61929866 

3. 27103405 

3. 87990075 
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Appendix  C 

DERIVATION  OF  Pp  FOR  EA  PROCESSOR  WITH  M  =  2  AND  ARBITRARY  N 


form 


For  the  EA  processor  with  M  =  2  and  arbitrary  N,  GUI  test  (24)  takes  the 


(C-l) 


here,  as  in  (B-l),  we  have  let  B^  =  |  and  Bj  =  jzj|.  The  threshold  y  need 
never  deviate  outside  the  range  [0,2];  the  lower  limit  on  the  range  follows 
because  the  left  side  of  (C-l)  can  never  be  negative.  The^ upper  limit  holds 
true  since  the  left  side  of  (C-l)  can  never  be  larger  than  2  E  B^  (by  Schwarz's 
inequality). 


For  signal-absent,  the  PDF  of  all  envelope  samples  in  (C-l)  is  Rayleigh: 

p(B)  =  B  exp(-B2/2),  B  >  0  .  (C-2) 

Since  a  common  absolute  scaling  of  envelopes  obviously  does  not  affect  the 
decision  made  in  (C-l),  a  convenient  scaling  has  been  selected  in  (C-2),  with 
no  loss  of  generality.  This  is  a  manifestation  of  the  CFAR  capability  of  test 
(C-l). 


Let  us  assume  temporarily  that  N  >  1  and  let 

N 

1  vS  ~2 

v  =  -  >  B.  . 

2  j 

j=l 

Then  the  PDF  of  random  variable  v  is 


(C-3) 


N-l  ,  t 

p(v)  =  X _ V  >  0 

'  (N  -  1) !  ’ 


(C-4) 


Now  let 


1/2  1/2 

B  =  (2r)  cos  6 ,  B  =  (2r)  sin  d  ,  0<r<oo,  0  <  8  <  r/2  ,  (C-5) 

X  L* 


C-l 
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in  (C-l),  and  employ  (C— 3) ;  there  follows  for  test  (C-l) 

sin(20)  >  7  —  +  7  -  1  . 

The  joint  PDF  for  random  variables  r  and  6  is  obtained  according  to 
p(r,0)  =  p(B  =  V2r  cos  d,  B  =  >/2r  sin  6)  x  |  Jacobian  | 

=  r  exp(-r)  sin(20),  0<r<oo,  0  <  9  <  tt/2  . 

Therefore  the  variables  r  and  6  are  independent  with  PDFs 

p(r)  =  r  exp(-r),  0  <  r  <  »  , 

p(0)  =  sin(2  0)  ,  0  <  $  <  t/2  . 

Now  define  the  random  variable 

q  =  v/r  . 


(C-6) 


{C-7) 


(C-8) 


(C-9) 


The  PDF  for  q  (which  is  independent  of  6)  is  given  by 
P(q)  =  J  dr  |r|  pr(r)  pv(qr) 

"~oo 


S' 


dr  r  r  exp(-r) 


(qr)N  1  exp(-qr) 
(N  -  1) ! 


(C-10) 


N  - 1 


=  N  (N  +  1) 


,N  +  2 


,  q  >  0  . 


(1  +  q) 

Furthermore,  test  (C-6)  can  now  be  expressed  as 

sin(20)  <  yq  +  y  -  1  . 


(C-ll) 


At  this  point  it  is  convenient  to  define  random  variable  u  =  sin(20). 
Then  the  probability  of  false  alarm  is  given  by 


P  =  Prob  (u  >yq  +  7  -  1) 

F 

/•OO  /•<*> 

=  J  dq  pq(q)  J  du  p^u)  . 
“°°  yq+7-1 


(C-12) 
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Now  the  inner  integral  in  (C-12)  can  be  evaluated  as  follows: 


r 

•'u 


du  Pu(u)  =  Prob  (u  >  U)  =  Prob  (sin (26)  >  U) 


/ 


=  2 


-ir/4 

l  d0  sin(20)  = 


—  arc  sin(U) 


1  , 


U  <  0 


2  1/2  (  (C-13) 

(1  -  U  )  ,  0  <  U  <  1 1 


0  , 


1  <  U 


Consequently  (C-12)  becomes 


1-7  2-7 

P^,  =  \  dq  p  (q)  1  +  j  7  dq  p  (q)  [l  -  (7q  +  7  -  l)2 

q  1-7  q 


I 


“|  l  /2 

J  • 


(C-14) 


At  this  point,  we  must  separate  the  derivation  into  two  parts:  one  for 
0  <  7  <  1  and  another  for  1  <  7  <  2 .  We  consider  the  latter  case  first. 
For  1  <  7  <  2  ,  (C-14)  becomes,  by  use  of  (C-10), 


2-7 

f  7  ,  q 

PF=J  dq  N  (N  +  1) - - 


N-  1 


'0 


(i  +  q) 


N  +  2 


[1  -  (yq  +  7  -  l)2] 


1/2 


2-7 

V 


=  N  (N  +  1)7  \  7  dq  Q 


N-  1 


(i  +  q) 


N  +  2 


(1  +q) 


1/2  ^2-7 


-  q 


1/2 


2-7 


Next  let  q  =  — - —  x  in  (C-15);  then  (reference  7,  15.  3. 1,  15. 1.  8) 


7 

-N-— 

^2  -7V’  2  r1  N-l  ~o  l  2-7  \  2 

P^,  =  N  (N  +  1)  7  ( - )  J  dx  x  (1  -  x)^  (1  + -  x] 


7 


1 

N+o  A 


7 


(C-15) 


=N(N  +  1)7(^l^lMF(N+|,N;N+|;.-i)  <c-. 


r(N  + 1) 

r(|)r(N  +  2)  i  n+1 

- jT-  (27)  2  (1  -  I)  2,  117  <2. 

r(N  +  2/ 
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Equation  (C- 16)  is  a  general  relation  for  Pp  for  y  in  the  range  [1,2].  Although 
(C-16)was  derived  for  N>1,  it  maybe  shown,  by  a  separate  derivation,  to  apply 
to  N  =  0  as  well.  (For  N  =  0,  the  allowed  range  for  y  does  not  include  (0, 1), 
but  is  restricted  to  [l,  2].) 


For  the  complementary  range  of  y ,  namely  [0,l],  we  return  to  (C-14) 
and  employ  (C-10)  to  obtain 


where 


1-7 


[1  =  j  7  dq  N 


PF=I1  +  I2’ 


N  -  1 


(N+  !)■ 


(i  +  q) 


N  +  2  ' 


r—  N-l  r  -1/2 

I  =  7  dq  N(N  +  1)  — *■■■  ;■-  [1  -  (yq  +  7  -  1)  J  * 

J  i-7  (i+q) 

y 

1  -  y 

In  integral  1^  ,  let  q= — — — x;  then  (reference  7,  15.3.1,  15.3.3) 


l1  =  N(N  + 


'KH1)  i0dxxN'1  (i+ifLxj 

-  (N  +  1)  (^j~)  F  +  2,  N;  N  +  1;  -  ^-) 


-N-2 


VN  .  -  v-N  - 1  -  v 

mn  +  1,(^2)  (1  +  L_1)  f(-i,  1;n  +  1;-^Z) 


(1  -  7)N  (N7  +  1)  . 


2/7 


In  integral  I2  ,  let  1  +  q  I  then  (reference  7,  15.3. 1) 

1 


i2  =  n(n*1>t(1-|)  I0  *  x2  (l  -  rr;  4 

,N-1 


N-l 


7  /  7\  /  3  5  7  \ 


(C-17) 


(C-18) 


(C-19) 


(C-20) 
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Combining  the  above  results  we  have 

PF  =  (l-T)N(l  +  N7)+N(N+l)i-(l-|)N  P(l-N,|if:  T^)’ 

(C-21) 

0  <  7  <  1  . 

The  final  results  for  Pf  are  given  by  (C-21)  and  (C-16)  for  the  complete 
range  [0,  2]  of  7.  In  table  C-l  the  required  values  of  threshold  y  are  presented. 


Table  C-l.  Threshold  Values  y 
for  EA  Processor  with  M  =  2 


N 

PF 

r 

0 

.01 

1.  999  949  998  75 

0 

.  001 

1.  999  999  500 

1 

.01 

1.  951  327  448 

1 

.  001 

1.989  581  459 

2 

.  01 

1.  796  631  253 

2 

.001 

1.  920  106  594 

3 

.01 

1.618  190  584 

3 

.  001 

1.805  310  077 

4' 

.  01 

1.454  698  805 

4 

.  001 

1.678  251  895 

6 

.  01 

1.  194  404  233 

6 

.  001 

1.442  868  746 

10 

.01 

.  867  343  947 

10 

.  001 

1.  100  639  981 

20 

.  01 

.  509  052  420 

20 

.  001 

.  676  673  936 
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Appendix  D 

SIMULATION  METHOD  AND  RESULTS  FOR  EA  PROCESSOR 


A  total  of  T  trials  is  conducted  for  the  random  variable  associated  with  the 
EA  processor  (see  (24)  and  (B-l)), 

/  M 

£•>) 

'  k=l  / 


M  N 

k=l  j=l 


(D-l) 


by  taking  G  groups,  each  of  the  size  S ,  where  T  =  GS .  In  the  g-th  group,  an 
estimate  qg  of  the  particular  quantile  of  interest  is  made,  based  upon  S 
samples.  The  sample  variance  of  the  G  quantile  estimates  is  then 


(D-2) 


Now  according  to  Cramer,  references,  sect.  28.5,  the  standard  devia¬ 
tion  of  a  quantile  estimate  based  upon  a  sample  of  size  S  is  proportional  to 
for  large  S.  We  estimate  the  proportionality  constant  K  by  setting 


(D-3) 


(D-4) 


The  estimate  of  the  quantile  from  the  total  of  T  trials  is  obtained  by 
combining  all  group  members  and  making  the  standard  selection  of  the  quantile. 
Its  standard  deviation  is  then  approximately 


(D-5) 
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Notice  that  we  have  not  used  a  Gaussian  assumption  for  the  quantile  estimates, 
but  only  the  inverse  square-root  decay  of  the  standard  deviation  of  the  quantile 
estimate  with  sample  size. 


The  random  variable  associated  with  the  EA  processor  in  (D-l)  is  com¬ 
pared  with  a  threshold  (scale  factor)  y  for  decisions  about  signal  presence  or 
absence.  We  let  y  be  the  estimated  threshold.  Table  D-l  presents  y  and  its 
estimated  standard  deviation.  (The  thresholds  for  M  =  1,  N  >  1  are  derivable 
from  the  threshold  table  for  the  UA  processor;  accordingly  no  results  for  this 
particular  case  are  presented  here. )  A  computer  program  for  the  estimation 
of  the  threshold  is  given  on  p.  D-4. 


As  N-»»  for  the  EA  processor,  GLR  test  (24)  tends  to 


£N 


k=l 


<  Threshold  . 


(D-6) 


The  characteristic  function  of  the  summation  random  variable  in  (D-6)  is 


dx  exp(i£x)  x  exp 


(D-7) 


This  quantity  can  be  evaluated  for  many  values  of  real  £  (for  any  M  and 
{dk})  via  M  fast  Fourier  transforms.  (For  equal  signal  amplitudes,  only  one 
FFT  is  needed. )  Then  Pd  is  available  via  the  techniques  in  reference  9. 
However,  this  technique  has  not  been  pursued  in  this  report;  hence  the  results 
for  N  =  °o  are  not  presented  in  figures  37-53  for  the  EA  processor. 
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Table  D-l.  Estimated  Thresholds  and  Standard  Deviations  for  the  EA  Processor 


P 

i-H 

O 

II 

Eli 

PF 

=  .001 

M 

N 

A 

Estimated  Standard 

A 

-y 

Estimated  Standard 

7 

Deviation  of  y 

1 

Deviation  of  y 

2 

0 

1.  99995"] 

-- 

1. 9999995' 

— 

1 

1.  95133 

— 

1. 98958 

— 

2 

1.  79663 

— 

1. 92011 

— 

3 

1.61819 

_ 

1. 80531 

_ 

>exact 

>exact 

* 

4 

1.45470 

— 

1. 67825 

— 

6 

1.  19440 

-- 

1.44287 

— 

10 

. 86734 

— 

1. 10064 

— 

20 

.  50905J 

— 

0.67667  w 

— 

3 

0 

2. 99414 

. 00036 

2. 99935 

. 00018 

1 

2. 8887 

.0029 

2.9660 

.  0032 

2 

2.6856 

.  0049 

2.8564 

.  0072 

3 

2.4600 

.  0072 

2. 6881 

.0093 

4 

2. 2574 

.0088 

2.552 

.  013 

6 

1.  9061 

.  0090 

2.245 

.  015 

10 

1.4452 

.  0100 

1.  781 

.  020 

20 

.  8874 

.0048 

1.  118 

.  016 

4 

0 

3. 9657 

.  0012 

3. 99152 

. 00070 

1 

3. 8052 

.  0045 

3. 9249 

.  0034 

2 

3.5351 

.0064 

3. 7849 

.  0088 

3 

3.3015 

.  0080 

3.  619 

.  012 

4 

3.0730 

.  0091 

3.430 

.  016 

6 

2.6587 

.  0124 

3.039 

.020 

10 

2.  0394 

.  0095 

2.454 

.  020 

20 

1.2996 

.  0059 

1.  608 

.018 

5 

0 

4.  9150 

.  0019 

4. 9756 

.  0018 

1 

4.6969 

.0041 

4. 8757 

.  0065 

2 

4.4389 

.  0069 

4. 6854 

.  0099 

3 

4.  1575 

.  0080 

4.497 

.  013 

4 

3. 8707 

.  0095 

4.272 

.023 

6 

3.4358 

.  0112 

3.  885 

.025 

10 

2.7378 

.  0106 

3.  151 

.  030 

20 

1.7744 

.  0101 

2.  180 

.  026 

10 

0 

9.4364 

.  0085 

9.685 

.  012 

1 

9.069 

.  010 

9.413 

.019 

2 

8.732 

.011 

9.  176 

.  019 

3 

8.404 

.  013 

8.  898 

.  026 

4 

8.\031\ 

.  015 

8.  541 

.  026 

6 

7.390 

.  018 

8.  006 

.031 

10 

6.319 

.  024 

6.  973 

.  036 

20 

4.514 

.  017 

5.  165 

.  042 
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PROGRAM  FOR  DETERMINING  QUANTILES 


PARAMETER  GROUPS=40,  MEMBRS=T000, TRIALS=GROUPS*MEMBRS 
INTEGER  OHM,OTM,OHT,OTT 
LOGICAL  REALS 

DIMENSION  R (TRIALS) ,  TH(GROUPS) ,  TT(GROUPS) 

DEFINE  IRAND=I*K+((1-SIGN(1,  I*K))/2)*34359738367 
DEFINE  RAND=FLOAT(I)/34359738367. 

M=10 

REALS=.  FALSE. 

K=5**15 

1=5281 

OHM=MEMBRS/lOO 
OTM=MEMBRS/1000 
OHT=TRIALS/100 
OTT  =TRIALS/ 1000 
DO  1  IA=0,  7 
N=IA 

IF  (IA.  EQ.  5)  N=6 

IF  (IA.  EQ.  6)  N=10 

IF  (IA.  EQ.  7)  N=20 

DO  21  ffil=l,  GROUPS 

ISTART=1+(IB1-1)*MEMBRS 

ISTOP=IBl*MEMBRS 

DO  2  IB2=1,  MEMBRS 

SL=0. 

SS=0. 

ST=0. 

DO  3  IC=1,  M 
I=IRAND 

RSQ=-LOG(l.  -RAND) 

SL=SL+SQRT  (RSQ) 

3  SS=SS+RSQ 

IF(N.  EQ.  0)  GO  TO  2 
DO  4  IC=1,  N 
I=IRAND 

RSQ=-LOG(l.  -RAND) 

4  ST=ST+RSQ 

2  R(ISTART-1+IB2)=SL*SL/ (SS+ST) 

CALL  SORT  (R,  1ST  ART ,  ISTOP,  REALS) 

TH(IB1)=.  5*(R(ISTOP-OHM)+R(ISTOP-OHM+l)) 
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21  TT(IB1)=.  5*(R(IST0P-0TM)+R(IST0P-0TM+1)) 

AVH=0. 

SDH=0. 

AVT=0. 

SDT=0. 

DO  22  IB1=1,  GROUPS 
A  VH=A  VH+T  H  (IB  1) 

AVT=AVT+TT  (IB1) 

SDH=SDH+TH(IB1)**2 

22  SDT=SDT+TT(IB1)**2 

SDH=SQRT((SDH-AVH**2/GROUPS)/((GROUPS-l)*GROUPS)) 
SDT=SQRT((SDT-AVT**2/GROUPS)/ ((GROUPS- l)*GROUPS)) 
CALL  SORT(R,  1,  TRIALS,  REALS) 

INTH=TRIALS-OHT 
INTT=TRIALS-OTT 
QH=.  5*(R(INTH)+R(INTH+1)) 

QT=.  5*(R(INTT)+R(INTT+1)) 

PRINT  23,  M,N,QH, SDH,QT,SDT 
PUNCH  23,  M,  N,  QH,  SDH,  QT,  SDT 

23  FORMAT  (21 1 0 ,  4E 1 5 .  9) 

1  CONTINUE 

END 
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